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ABSTRACT 

We carefully consider observational and theoretical constraints on the global properties of 
secondary stars in cataclysmic variable stars (CVs). We then use these constraints to construct 
and test a complete, semi-empirical donor sequence for CVs with orbital periods Porb ^ 6 hrs. 
All key physical and photometric parameters of CV secondaries (along with their spectral 
types) are given as a function of Porb along this sequence. This provides a benchmark for 
observational and theoretical studies of CV donors and evolution. 

The main observational basis for our donor sequence is an empirical mass-radius rela- 
tionship for CV secondaries. Patterson and co-workers have recently shown that this can be 
derived from superhumping and/or eclipsing CVs. We independently revisit all of the key 
steps in this derivation, including the calibration of the period excess-mass ratio relation for 
superhumpers and the use of a single representative primary mass for most CVs. We also 
present an optimal techique for estimating the parameters of the mass-radius relation that si- 
multaneously ensures consistency with the observed locations of the period gap and the period 
minimum. We present new determinations of these periods, finding Pgap,+ — 3.18 ± 0.04 hrs 
(upper edge), Pgap.- = 2.15 ± 0.03 hrs (lower edge) and Pmin = 76.2 ± 1.0 min (period 
minimum). 

We test the donor sequence by comparing observed and predicted spectral types (SpTs) 
as a function of orbital period. To this end, we update the SpT compilation of Beuermann and 
co-workers and show explicitly that CV donors have later SpTs than main sequence (MS) 
stars at all orbital periods. This extends the conclusion of the earlier study to the short-period 
regime (Porb < 3 hrs). We then compare our donor sequence to the CV data, and find that it 
does an excellent job of matching the observed SpTs. Thus the empirical mass-radius relation 
yields just the right amount of radius expansion to account for the later-than-MS spectral 
types of CV donors. There is remarkably little intrinsic scatter in both the mass-radius and 
SpT — Porb relations, which confirms that most CVs follow a unique evolution track. 

The donor sequence exhibits a fairly sharp drop in temperature, luminosity, and opti- 
cal/infrared flux well before the minimum period. This may help to explain why the detection 
of brown dwarf secondaries in CVs has proven to be extremely difficult. 

We finally apply the donor sequence to the problem of distance estimation. Based on 
a sample of 22 CVs with trigonometric parallaxes and reliable 2MASS data, we show that 
the donor sequence correctly traces the envelope of the observed M.jhk ~ Porb distribution. 
Thus robust lower limits on distances can be obtained from single-epoch infrared observa- 
tions. However, for our sample, these limits are typically about a factor of two below the true 
distances. 

Key words: accretion, accretion disks - stars: novae, cataclysmic variables - stars: distances 
- stars: fundamental parameters. 



1 INTRODUCTION 

Cataclysmic variable stars (CVs) are compact, interacting binary 
systems in which a white dwarf primary accretes from a low-mass, 
roughly main-sequence donor star. The mass transfer and secular 
evolution of CVs is driven by angular momentum losses. In systems 
with long orbital periods (Port <^ 3 hrs), the dominant angular mo- 



mentum loss mechanism is thought to be magnetic braking (MB) 
due to a stellar wind from the donor star. In the canonical "dis- 
rupted magnetic braking" evolution scenario for CVs, MB stops 
when the secondary becomes fully convective, at Port — 3 hrs. 
At this point, the donor detaches from the Roche lobe, and grav- 
itational radiation (GR) becomes the only remaining angular mo- 
mentum loss mechanism. The GR-driven shrinkage of the orbit ul- 
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timately brings the secondary back into contact at Port — 2 hrs, 
at which point mass transfer resumes. The motivation for this sce- 
nario is the dearth of mass-transferring CVs in the so-called period 
gap between Porb 2 hrs and Porb — 3 hrs. 

As a CV evolves, its donor star is continuously losing mass. As 
long as the mass-loss time scale (tj^j^ ~ M2/M2) is much longer 

than the donor's thermal time scale (rth ~ 'r^l^ ^' secondary 
is able to maintain thermal equilibrium and should closely follow 
the standard main sequence track defined by single stars. Here and 
throughout, we use M2, R2 and L2 to denote a donor's mass, radius 
and luminosity, and M2 to denote the rate at which it is losing mass 
to the primary. If the condition Tj^^ » Tm^^ is not met, the donor 
will be driven out of thermal equilibrium and become oversized 
compared to an isolated main sequence (MS) star of identical mass. 

So how do thermal and mass-loss time scales compare for the 
donor stars in CVs? Or, to put it another way, are the donors main 
sequence stars? This question has been addressed twice in recent 
years, in very different ways. First, Beuermann et al. (1998; here- 
after B98) showed that, at least for periods Porb ^ 3 hrs (i.e. above 
the period gap), the spectral types (SpTs) of CV donor stars are sig- 
nificantly later than those of isolated MS stars. ^ At the very longest 
periods (Port ^5 — 6 hrs), this is probably due to the donors being 
somewhat evolved (B98; also see Podsiadlowski, Han & Rappaport 
2003). However, for all other systems, the late SpTs of the donors 
are simply a sign of their losing battle to maintain thermal equilib- 
ritun. 

Second, Patterson et al. (2005; hereafter P05) produced an em- 
pirical mass-radius sequence for CV donors based on a sample of 
masses and radii derived mainly from superhumping CVs. One of 
their key results was that C V donors are indeed oversized compared 
to MS stars, for all periods Porb ^ 6 hrs. They also detected a clear 
discontinuity in R2 for systems with similar M2 above and below 
the period gap. This is just what is expected in the disrupted mag- 
netic braking picture: the MB-driven systems just above the gap 
should be losing mass faster than than the GR-driven systems just 
below. They should therefore be further out of thermal equilibrium. 

In the present paper, we will build on these important stud- 
ies by constructing a complete, semi-empirical donor sequence for 
CVs. More specifically, our goal is to derive a benchmark evolution 
track for CV secondaries that incorporates all of the best existing 
observational and theoretical constraints on the global donor prop- 
erties. ^ Along the way, we will update both of the earUer studies 
and show that they are mutually consistent: the mass-radius rela- 
tion obtained by P05 is just what is needed to account for the late 
SpTs observed by B98. 

Our donor sequence should be useful for many practical ap- 
plications. Perhaps most fundamentally, it provides a benchmark 
for what we mean by a "normal" CV and can thus be used to test 
CV evolution scenarios that predict the properties of CV donors. 
The sequence can also be used to simplify and improve on "Bai- 
ley's method" for estimating the distances to CVs (Bailey 1981). 
The new version of the method requires only knowledge of the or- 



^ For the record, earlier studies along these lines were carried out by 
Echevam'a (1983), Patterson (1984), Friend et al (199()ab) and Smith & 
Dhillon (1998). 

^ Here and below, we use the term "donor sequence" to describe the de- 
pendence of the secondary star parameters on the orbital period of a CV. 
We call a sequence "complete" if aU important physical {Porb> ^2, R2, 
T^ff, SpT) and photometric (absolute magnitudes in UBVRUHKLM) 
donor properties are fully specified along the entire evolutionary track. 



bital period and an infrared magnitude measurement. Conversely, 
for systems with known distances, the sequence can be used to es- 
timate the donor contribution to the system's flux in any desired 
bandpass. Finally, the present work is also a stepping stone towards 
another goal: the construction of new, semi-empirical mass-transfer 
rate and angular momentum loss laws for CVs. This should be pos- 
sible, since the degree of departure of a donor from the MS track is 
a direct measure of the mass loss from it and hence of the angular 
momentum loss that drives M2. 



2 THE EMPIRICAL MASS-RADIUS RELATION FOR CV 
DONOR STARS 

Given the importance of the M2 — R2 relation to the present work, 
we begin by revisiting and updating P05's M2-R2 relation for CV 
donor stars. We still use the same fundamental data as P05, but 
carry out a fuUy independent analysis to derive our ovm mass- 
radius relation from the data. In the following sections, we will 
briefly discuss each of the key steps in the derivation (and also in- 
troduce some new ideas of our own). However, we start with an 
outline of the overall framework of the method, i.e. the way in 
which a mass-radius relation for CV secondaries can be derived 
from (mainly) observations of superhumps. 

2.1 Donor Masses and Radii from Superliumping CVs 

Superhumps are a manifestation of a donor-induced accretion disk 
eccentricity that is observed primarily in erupting dwarf novae, 
but also in some non-magnetic nova-like CVs and even in a few 
LMXBs. Once established, the eccentricity precesses on a time 
scale that is much longer than the orbital period. The superhump 
signal is then observed at the beat period between the orbital and 
precession periods. The superhump period, Psh, is therefore t5'p- 
ically a few percent longer than the orbital period, and the super- 
himip excess, e, is defined as 



Both theory and observation agree that e is a function of the 
mass ratio, q — M2/M1. This e ~ q relation can be calibrated 
by considering eclipsing superhumpers for which an independent 
mass ratio estimate is available. Given such a calibration, we can 
estimate q for any superhumper with measured e. 

Two additional steps are needed to turn an estimate of the mass 
ratio into an estimate of Ad2 and i?2. First, in order to obtain M2 
from q, we clearly first need an estimate of Mi. For a few sys- 
tems (mostly eclipsers), this can be measured directly. However, 
for most other systems, a representative value has to be assumed. 
For example, P05 used Mi = O.TSM© for all systems without a 
direct estimate, based on several estimates of the mean WD mass in 
CVs in the literature. With this assumption, the mass of the donor 
is then simply given by 

M2=qMi. (2) 

The second additional step is to make use of the fact that the 
secondary is filling its Roche lobe. The donor must therefore obey 
the well-known period-density relation for Roche-lobe filling ob- 
jects (Warner 1995, Equation 2.3b) 

< p2 >= 107 P-If, g cm^ (3) 
where Porb,h is the orbital period in units of hours and 
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< P2 > = 



3M2 



(4) 



Equation 3 is accurate to better than 3% over the interval 0.01 < 
g < 1. Since M2 and Port are known at this point, the period- 
density relation can be recast to yield an estimate of the donor ra- 
dius 



R2/Re = 0.2361 P^^l^ q'^^ (M^/Mq)'^' 
or, equivalently, 

R2/RQ = 0.2361 P^ll,^ (Afa/Mo)^/^ 



(5) 



(6) 



At this point, both M2 and R2 have been estimated. A fit to the 
resulting set of M2 and R2 pairs (supplemented with additional 

data points derived from from eclipsing systems) can then be used 
to determine the functional form of the mass-radius relation for CV 
donors. 

In the following sections, we will take a closer look at all of the 
key steps we have outlined, and also add one final step of our own. 
More specifically, we will consider (i) the calibration of the e — q 
relation for superhumpers; (ii) the assumption of a single Mi value 
for most CVs; (iii) the external constraints that the final mass-radius 
relation should reproduce the observed locations of the period gap 
and the minimum period; (iv) the derivation of an optimal fit to 
the M2-R2 pairs, allowing for correlated errors, intrinsic dispersion 
and external constraints. 



2.2 Calibrating the e — q Relation 

Table 7 in P05 provides a list of calibrators for the e — g relation. 
This contains 10 superhumping and eclipsing CVs with indepen- 
dent mass ratio constraints, one superhumping CV with a large su- 
perhump excess and an assumed upper limit on q (BB Dor), and 
also one superhumping and eclipsing LMXB with a very low mass 
ratio (KV UMa). In devising our own, independent calibration of 
the e — q relation, we used the same set of calibrators, but analysed 
them independently. Details are given in Appendix A; the results 
are shown in Figure 1. Our preferred fit to the data is given by 



g(e) = (0.114 ± 0.005) + (3.97 ± 0.41) x (e - 0.025) 



(7) 



The errors here are l-a for both parameters jointly, and the shift 
applied to e ensures that the fit parameters (and their errors) are 
reasonably uncorrelated. The fit is shown in Figure 1 and achieves 
a statistically acceptable xt = 103 without the need to add any 
intrinsic dispersion in excess of the statistical errors on the e and q 
estimates. Thus any intrinsic scatter around the calibrating relation 
must be small compared to the statistical errors on the data points. 
Figure 1 also allows a direct comparison of our fit to the data against 
P05's, as well as against two recent theoretically-motivated calibra- 
tions of the e — q relationship (Goodchild & Ogilvie 2006; Pearson 
2006). All calibrations agree quite well, except at the highest mass 
ratios, where data are sparse. 

One final point worth noting is that the statistical errors on 
the fit parameters (and indeed the uncertainty regarding the func- 
tional form of the calibration itself) translate into systematic errors 
on the resulting masses and radii. The l-a error band arising from 
the statistical uncertainties on our fit parameters is indicated by the 
dashed region in Figure 1. Formally, this is less than 10% even out 
to Q ~ 0.4, but the fit is very poorly constrained beyond q ~ 0.3. 
The e—q relation could thus change shape in this regime. The statis- 
tical error on a mass ratio obtained via Equation 7 can be estimated 
by folding the error on e through the calibrating relation. 




Figure 1. Calibration of the q—c relation for superhumping CVs. Top panel: 
Constraints on the slope and intercept of the preferred linear calibration. The 
cross marks the best-fitting parameter combination. The eUipses correspond 
to la, 2a and 3a contours on both parameters jointly. Bottom panel: Mass 
ratio vs period excess for the calibrating stars. The thick solid line is the 
best-fitting linear calibration, with the shaded area marking the statistical 1- 
a range. The thick long-dashed line shows P05's calibration, the thin dotted 
line is Pearson's (2006) calibration, and the thin short-dashed line is the 
calibration of Goodchild & Ogilvie (2006). The data point shown as a thin 
open symbol shows the upper limit on the mass ratio of BB Dor proposed 
by P05. 



2.3 The Assumption of Constant Primary Mass 

The next key step in the derivation of masses and radii from su- 
perhumpers is the assumption of a single primary mass for most 
systems in the sample. The worry here is not so much that the as- 
sumed mass might differ from the true mean WD mass for CVs. 
This would simply shift all (log M2, log R2) pairs along a line of 
slope 1/3, but would not alter the shape of the donor mass-radius 
relation. Instead, the main concern is that Mi might exhibit system- 
atic trends within the observed CV population, most importantly 
with orbital period. Such trends could affect the shape of the mass- 
radius relation. 

In order to test whether this is a problem, we can check if there 
is any dependence of M\ on Porb among CVs with reliable WD 
mass estimates (i.e. eclipsing systems). P05 actually tabulate all 
available Mi estimates derived from eclipsing CVs. In the bottom 
panel of Figure 2, we plot these estimates as a function of Port for 
all systems with Porb < 6 hrs. We exclude systems with longer 
periods, since they have evolved donors and thus follow a different 
evolution track (see Section 3.3 below). 

Figure 2 does not provide evidence for evolution of Mi with 
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Figure 2. Constraints on white dwarf mass evolution with orbital period. 
Top panel: Results of a linear fit to Mw D vs Porb for eclipsing systems 
with weU-determined white dwarf mass. The cross marks the best-fitting 
parameter combination. The ellipses correspond to la, 2a and 3a contours 
on both parameters jointly. Bottom panel: My/ jj vs Porb for eclipsers. The 
thick solid line shows the optimally estimate mean white dwarf mass for all 
systems: My^u = 0.75 ± 0.05Mq The thin dashed line shows the best 
linear fit to the data. 



Porb- A full discussion of the statistical evidence for this assertion 
is given in Appendix B, along with a discussion of how this state- 
ment can be reconciled with earlier work that seemed to show such 
evolution. Here, we simply note that a linear fit to the data yields 
a slope consistent with zero. The best-fitting constant Mi yields a 
mean WD mass of < Mi 0.75 ± O.OSMq and requires an in- 
trinsic scatter of atnt = O.I6M0. Thus we confirm P05's assertion 
that, in the absence of other information. Mi can be estimated as 
Ml = O.75M0 with about 20% uncertainty. 

We finally note that the global intrinsic dispersion - ctmi = 
0.16 Mq, i.e. 21% - is the statistical error associated with taking 
Ml = O.75M0 for any particular system. By contrast, the formal 
error on the mean WD mass (O.OSMq) translates into a systematic 
error on donor masses and radii (since a shift in the assvuned Mi 
affects all affected data points in the same way). 

Having dealt with the e — q calibration and the assumption 
of constant Mi, we are ready to calculate an updated set of M2 
and R2 values (with statistical errors) for all superhumpers. We 
list these in Table 1. Our estimates are not identical to those de- 
rived by P05, since we have used a different e — q calibration and 
have explicitly accounted for the intrinsic dispersion in Mi when 
estimating the statistical errors. However, as expected, the overall 
pattern presented by the data is much the same. The main change in 



the data values is a slight shift towards higher masses, which arises 
because our g-estimates are generally a little higher than P05's for 
fixed e (see Figure 1). 



2.4 External Constraints: The Location of the Period Gap 
and the Minimum Period 

In principle, wc could now simply fit the M2 — R2 pairs in Table 1 
(supplemented with data for eclipsing systems). However, there are 
actually additional empirical constraints that can (and should) be 
imposed on the mass-radius relation. These constraint come from 
the observed locations of the period gap and of the period mini- 
mum. 

Let us first consider the period gap. The bottom panel in Fig- 
ure 3 shows the M2 — R2 estimates from Table 1, along with sim- 
ilar data for eclipsing systems, taken from Table 8 in P05. There 
is a clear discontinuity in donor radii at M2 — O.2M0, with 
donors in long-period systems being larger than those in short- 
period systems. As discussed in more detail by P05, the transition 
is quite sharp and broadly consistent with the standard CV evo- 
lution scenario. Briefly, we should expect systems just above and 
below the period gap to have identical donor masses, since CVs 
evolve through the period gap as detached binaries, with no sig- 
nificant mass loss from the secondary. However, their radii must 
differ, since both short- and long-period donors must still obey the 
period-density relation for Roche-lobe- filling stars (Equation 3). 
More specifically, if we denote the upper and lower edges of the 
period gap as Pgap,±, the ratio of donor radii at the gap edges must 
satisfy 



R2, 



R2,- 



2/3 



(8) 



Physically, donors above the gap are larger since they lose mass at 
a much higher rate than those below and are therefore forced more 
out of thermal equilibrium. 

If we accept the premise that CV donors evolve through the 
period gap as detached systems, the data in Figure 3 suggests that 
the critical donor mass at which mass transfer stops (and then 
restarts) is Mconv = 0.20 ± O.O2M0, where the estimated error is 
purely statistical. ^ In other words, Mconv is the donor mass at both 
edges of the period gap. Given empirical estimates for the location 
of these edges, we can therefore use the period-density relation to 
determine Ri.conv, i.e. the radii of donors at the gap edges. Thus 
the locations of the gap edges iPgap,±) fix the A/2 — R2 relations 
for long- and short-period systems at Mconv and Pgap,±- 

Figure 4 shows the orbital period distribution of CVs drawn 
from the Ritter & Kolb (2003) catalog (Edition 7.6) in both dif- 
ferential and cumulative form. The period gap is obvious in these 
distributions. We have carried out repeated measurements of the 
gap edges from this data with a variety of binning schemes. Based 



on these measurements, we estimate Pn, 



= 2.15±0.03 hrsand 



Pgap,+ = 3.18 ± 0.04 hrs. It is worth noting that the number of 
systems inside the gap increases towards longer periods. This may 
be due to low-metallicity systems, which are expected to invade the 
gap from above (Webbink & Wickramasinghe 2002). 

We can similarly demand that our mass-radius relationship 
should reproduce the observed minimum period, Pmin, of the CV 

^ We use the subscript com to denote this critical mass, since in the canon- 
ical picture it corresponds to the mass at which the donor becomes fully 
convective. 
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Figure 3. The mass-radius relation of CV donor stars. Bottom panel: Points shown are empirical mass and radius estimates for CV donors. Superhumpers 
are shown in black, eclipsers in red. Filled squares correspond to short-period CVs, filled circles to long-period systems, and crosses to likely period 
bouncers. The parallelograms illustrate the typical error on a single short-period or long-period CV. Open symbols were ignored in fits to the data 
since they correspond to systems in the period gap or long-period (probably evolved) systems. The solid lines show the optimal fit to the data in the 
period bouncer, short-period and long-period regions. The dotted fine is the mass-radius relation for main sequence stars taken from the 5 Gyr BCAH98 
isochrone. Top panels: Constraints on the power law exponents of the M2 — R2 relations in the three period/mass regimes. For each regime, we plot 

vs exponent and indicate the corresponding to l-a, 2-a and 3-a around the minimum with horizontal dashed lines. The strange shape of the 
curve for the period bouncers near the exponent i is real. It arises because the intrinsic power law relationship between M2 and R2 estimates has the 
same exponent. 
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Table 1. Donor masses and radii as estimated from superiiump periods. Orbital periods are in hours; masses and radii are in solar units. Note that mass and 
radius estimates are not independent, but correlated via Equation 6. 



System 


^orb 


Ma 


0"M2 


R2 




System 


^orh 


Ma 




R2 




DIUMa 


1.3094 


0.051 


0.011 


0.105 


0.008 


CYUMa 


1.6697 


0.119 


0.026 


0.164 


0.012 


V844 Her 


1.3114 


0.083 


0.018 


0.124 


0.009 


FO And 


1.7186 


0.115 


0.027 


0.165 


0.013 


LL And 


1.3212 


0.097 


0.023 


0.131 


0.010 


OUVir 


1.7450 


0.130 


0.029 


0.173 


0.013 


SDSS 0137-09 


1.3289 


0.085 


0.019 


0.125 


0.009 


VZPyx 


1.7597 


0.110 


0.024 


0.165 


0.012 


ASAS 0025+12 


1.3452 


0.072 


0.017 


0.120 


0.009 


CC Cnc 


1.7645 


0.156 


0.034 


0.186 


0.014 


AL Com 


1.3601 


0.047 


0.010 


0.104 


0.008 


HT Cas 


1.7676 


0.089 


0.009 


0.154 


0.005 


WZSge 


1.3606 


0.056 


0.009 


0.111 


0.006 


lYUMa 


1.7738 


0.093 


0.006 


0.157 


0.003 


RX 1839+26 


1.3606 


0.063 


0.015 


0.115 


0.009 


VWHyi 


1.7825 


0.110 


0.024 


0.166 


0.012 


PUCMa 


1.3606 


0.077 


0.018 


0.123 


0.009 


ZCha 


1.7880 


0.089 


0.003 


0.155 


0.001 


SWUMa 


1.3634 


0.084 


0.020 


0.127 


0.010 


QWSer 


1.7887 


0.110 


0.026 


0.167 


0.013 


HV Vir 


1.3697 


0.071 


0.015 


0.120 


0.009 


WXHyi 


1.7954 


0.114 


0.025 


0.169 


0.012 


MM Hya 


1.3822 


0.066 


0.014 


0.118 


0.009 


BKLyn 


1.7995 


0.154 


0.033 


0.187 


0.013 


WXCet 


1.3990 


0.070 


0.016 


0.122 


0.009 


RZLeo 


1.8250 


0.114 


0.026 


0.171 


0.013 


KVDra 


1.4102 


0.080 


0.018 


0.128 


0.010 


AW Gem 


1.8290 


0.137 


0.030 


0.182 


0.013 


T Leo 


1.4117 


0.081 


0.018 


0.129 


0.009 


SU UMa 


1.8324 


0.105 


0.023 


0.167 


0.012 


EGCnc 


1.4393 


0.031 


0.007 


0.095 


0.007 


SDSS 1730+62 


1.8372 


0.123 


0.027 


0.176 


0.013 


V1040 Cen 


1.4467 


0.103 


0.023 


0.142 


0.011 


HS Vu- 


1.8456 


0.153 


0.033 


0.190 


0.014 


RXVol 


1.4472 


0.064 


0.015 


0.121 


0.009 


V503 Cyg 


1.8648 


0.139 


0.031 


0.185 


0.014 


AQEri 


1.4626 


0.096 


0.021 


0.139 


0.010 


V359 Cen 


1.8696 


0.127 


0.030 


0.180 


0.014 


XZEri 


1.4678 


0.094 


0.005 


0.139 


0.003 


CUVel 


1.8840 


0.098 


0.024 


0.166 


0.013 


CP Pup 


1.4748 


0.083 


0.015 


0.133 


0.008 


NSV 9923 


1.8984 


0.134 


0.030 


0.185 


0.014 


V1159 0ri 


1.4923 


0.106 


0.023 


0.146 


0.010 


BRLup 


1.9080 


0.112 


0.027 


0.175 


0.014 


V2051 Oph 


1.4983 


0.095 


0.022 


0.141 


0.011 


VI 974 Cyg 


1.9502 


0.202 


0.033 


0.216 


0.012 


V436 Cen 


1.5000 


0.074 


0.018 


0.130 


0.011 


TUCrt 


1.9702 


0.129 


0.028 


0.188 


0.014 


BCUMa 


1.5026 


0.102 


0.022 


0.145 


0.010 


TYPsA 


2.0194 


0.135 


0.030 


0.194 


0.014 


HO Del 


1.5038 


0.093 


0.022 


0.141 


0.011 


KKTel 


2.0287 


0.121 


0.027 


0.187 


0.014 


EKTrA 


1.5091 


0.107 


0.024 


0.147 


0.011 


V452 Cas 


2.0304 


0.159 


0.035 


0.205 


0.015 


TVCrv 


1.5096 


0.108 


0.025 


0.148 


0.011 


DVUMa 


2.0604 


0.165 


0.013 


0.209 


0.006 


VY Aqr 


1.5142 


0.072 


0.016 


0.129 


0.010 


YZCnc 


2.0832 


0.176 


0.038 


0.216 


0.016 


OYCar 


1.5149 


0.065 


0.004 


0.125 


0.003 


GXCas 


2.1365 


0.145 


0.032 


0.206 


0.015 


RX 1131+43 


1.5194 


0.088 


0.019 


0.139 


0.010 


NY Ser 


2.3460 


0.197 


0.043 


0.242 


0.018 


ERUMa 


1.5278 


0.105 


0.023 


0.148 


0.011 


V348 Pup 


2.4442 


0.202 


0.045 


0.251 


0.019 


DMLyr 


1.5710 


0.095 


0.022 


0.145 


0.011 


V795 Her 


2.5982 


0.237 


0.051 


0.276 


0.020 


UVPer 


1.5574 


0.081 


0.019 


0.137 


0.010 


V592 Cas 


2.7614 


0.197 


0.042 


0.270 


0.019 


AKCnc 


1.5624 


0.121 


0.028 


0.157 


0.012 


TUMen 


2.8128 


0.225 


0.049 


0.286 


0.021 


AO Oct 


1.5737 


0.083 


0.021 


0.139 


0.012 


AH Men 


3.0530 


0.275 


0.059 


0.323 


0.023 


SXLMi 


1.6121 


0.114 


0.026 


0.158 


0.012 


DWUMa 


3.2786 


0.197 


0.010 


0.303 


0.005 


SSUMi 


1.6267 


0.118 


0.026 


0.160 


0.012 


TTAri 


3.3012 


0.263 


0.056 


0.335 


0.024 


KSUMa 


1.6310 


0.083 


0.020 


0.143 


0.011 


V603 Aql 


3.3144 


0.242 


0.044 


0.327 


0.020 


V1208 Tau 


1.6344 


0.122 


0.027 


0.163 


0.012 


PX And 


3.5124 


0.278 


0.060 


0.356 


0.025 


RZSge 


1.6387 


0.102 


0.023 


0.153 


0.012 


V533 Her 


3.5352 


0.225 


0.048 


0.333 


0.024 


TY Psc 


1.6399 


0.114 


0.025 


0.159 


0.012 


BE Dor 


3.5808 


0.291 


0.062 


0.366 


0.026 


IRGem 


1.6416 


0.116 


0.032 


0.160 


0.015 


BHLyn 


3.7380 


0.246 


0.053 


0.356 


0.026 


V699 Oph 


1.6536 


0.070 


0.017 


0.136 


0.011 


UUAqr 


3.9259 


0.197 


0.041 


0.342 


0.024 



population. Based again on the data in Figure 4, we estimate this to 
be Pmin — 76.2 lb 1.0 min. We can implement this constraint by 
truncating the M2 — R2 relation for short-period CVs at Mhounce, 
the donor mass where it reaches Pmin- If there are data points with 
lower masses, they need to be fit separately, subject to the constraint 
that the fit should meet the mass-radius relation for short-period 
CVs at Mbounce- Moreovcr, if Pmin is supposed to be a minimum 
period, then the fit to systems with M2 < Mhounce should yield an 
M2 — R2 relation that corresponds to increasing orbital period with 
decreasing M2. We will find below that that is indeed the case. 



the overall M2 — R2 relation? This question is not as trivial as 
it may seem, for two reasons. First, the external constraints (i.e. 
the measured values of Pmin, Pgap,± and Mconv) need to be im- 
posed self-consistently on fits to the M2 — R2 pairs. Second, the 
mass and radius estimates for any given superhumper (along with 
their errors) are completely correlated. This was already pointed 
out by P05 and is easy to see from Equation 6, which shows that 

1 /3 

R2 oc Mj ' . The masses and radii of eclipsers are similarly corre- 
lated, since only q and Mi are generally estimated directly from the 
light curve, with M2 and R2 being obtained indirectly in much the 



2.5 The Optimal M2 -R2 Relation for C V Donors 

Given our set of mass-radius pairs and the added external con- 
straints, how can we obtain an optimal analytical description of 
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2 3 4 

Orbital Period (hrs) 



Figure 4. Differential and cumulative orbital period distribution of CVs, 
based on data taken from Edition 7.6 of the Ritter & Kolb (2003) catalog. 
Our adopted values for the minimum period and the period gap edges are 
shown as vertical lines. The shaded regions around them indicate our esti- 
mate of the errors on these values. 



same way as for superhumpers. This correlation should be taken 
into account when fitting the data. 

In Appendix C, we show how the standard x^-statistic can be 
modified so as to explicitly account for correlated data points, exter- 
nal constraints and intrinsic dispersion. When carrying out our 
fits, we explicitly distinguish between three groups of CVs, each of 
which is fit independently with a power-law M2 — R2 relation: (i) 
non-evolved long-period CVs (3.18 hrs < Port < 6 hrs); (ii) "nor- 
mal" short-period CVs (Port < 2.15 hrs and M2 > Mtounce', 
(iii) period bouncers (A/2 < Altounce)- The value of Mtounce is 
obtained from the fit to short-period CV, as described in Section 2.4. 
In principle, this could be an iterative procedure, since the value of 
Mbounce determines which points are included in the fit. In prac- 
tice, the fit to the short-period systems is quite insensitive to the 
exact value, since the vast majority of points lie well above any 
plausible estimate for Mbounce- For the optimal fit described be- 
low, we find Mbounce = O.O626Af0. 

Remarkably, we find that the intrinsic radius scatter required 
by the }^ fits to the short- and long-period systems is only about 
2%-3%. For comparison, the typical statistical error is about 7%- 
8%. This implies that the intrinsic dispersion is probably rather 
poorly constrained by the data, but also that it cannot be more than a 
few percent. Thus there really is a unique evolution track that is fol- 
lowed by almost all CVs. When other tracers are used to study CV 
evolution (e.g. luminosity or period changes), the existence of this 
track is often masked by variability on time scales shorter than the 
binary evolution time scale (e.g. Biining & Ritter 2004). However, 



* We actually calculated our own R2 estimates and eiTors for echpsers 
from the M\, M2 and Porb values listed in Table 8 of P05; this ensures 
consistency, i.e. all of the data satisfy the same period-density relation. As 
expected, our numbers agreed very well with those in P05's Table 8. 



donor radii are hardly affected by these short-time-scale effects and 
thus faithfully trace the long-term evolution of M2 

Our final mass-radius relation for all three types of CV donors 
is given by: 



Rq 



(0.110 ± 0.005) X (^ T^^^ y 



, +0.05 
-0.10 



[0.230 ± 0.008) X 



(9) 



(0.299 ± 0.010) X 



where the three regimes are formally defined as (i) period bouncers: 

M2 < Mbounce; (ii) short-period CVs: Mbounce < M2 < Mconv 

andPori, < Pgap,-; (iii) long-period CVs: Mconv < M2 < M^^oi 



and P„. 



< Porb < Pevoi- For reference, let us also summarize 



the various quantities that have been assumed 



^^bounce 


= 0.063 ± 0.009 Mq 


^^conv 


= 0.20 ± 0.02 Mq 


Mevol 


~ 0.6- 0.7 Mq 


^ rmn 


= 76.2 ± 1.0 min 


P 


= 2.15 ± 0.03 hr 


Pgap,+ 


= 3.18 ± 0.04 hr 


Pevol 


~ 5 - 6 hr. 



(10) 



We note in passing that the best-fitting mass-radius relation 
in the short-period regime is slightly steeper than might be sug- 
gested by visual inspection of the data in Figure 3. This is because 
M2 and R2 are correlated. More specifically, the appearance of the 
power law index b in Equation C5 acts to push the fit away from 
indeces 6 — | and thus in this case to larger values. This can be 
verified by artificially increasing the assumed intrinsic dispersion, 
since this forces the fit towards a standard least-squares solution. 
As expected, the power law index then tends to a slightly shallower 
value of 0.60 in the short-period regime, and the donor mass at pe- 
riod minimum becomes Mbounce ~ O.O53A'/0. 

P05 did not derive a fit to the period bouncers, due to the 
sparseness of the donor data in the sub-stellar regime. We com- 
pletely agree with P05 that the fit is very poorly constrained for 
this group: only 5 data points fall into this class in our sample, and 
only one of these is clearly inconsistent with the A/2 — R2 relation 
for "ordinary" short-period systems. The A/2 — R2 relation in this 
regime is thus strongly affected by our estimate for Pmin and by 
the parameters inferred for short-period systems. These constraints 
ultimately decide which points are included in the period bouncer 
class and also uniquely fix R2 at Mbounce- We have nevertheless 
provided the best-fit parameters for this group, partly because we 
have no reason to distrust these constraints, and partly because we 
want our semi-empirical donor sequence to extend into the sub- 
stellar (period bounce) regime; we would therefore rather have an 
uncertain mass-radius relation than none. 

We also derive some confidence in our fit to the period bounc- 
ers from the fact that the best-fit mass-radius exponent in this 
regime turns out to be less than 1/3 (at somewhat better than 2- 
cr). This is a critical value, because the effective mass-radius index 
of the donor determines the sign of the period derivative of a CV. 
This can be seen explicitly by combining Paczyniski's (1971) ap- 
proximation for the radius of Roche-lobe-filling secondaries 



^ = 0.462 ( ^ 

a \l + q 



1/3 



(11) 



(which is valid for q ^ 0.8) with Kepler's third law to obtain 
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P'orb ^ 

Porb 2 Ma' ^ ' 

Here, M2 < is the mass transfer rate from the donor, and ^ is 
the effective mass-radius index of the donor along its evolutionary 
track. If the mass-radius relation along the track can be described 
as a power law, ^ is simply the power-law exponent. Equation 12 
shows that the orbital period decreases for ^ > 1/3, but increases 
for ^ < 1/3. With = 0.21+" "ij for our period bouncers, these 
systems appear to be evolving back towards larger periods - as they 
should, since Pmm is supposed to be a minimum period along our 
donor track. 



3 THE LATE SPECTRAL TYPES OF CV DONORS 

We now turn our attention to the spectral types of CV donors, and 
specifically to the SpT-Porb relationship. Ultimately, our goal will 
be to test if donor models based on the M2 — R2 relations we have 
just derived are capable of reproducing the empirical SpT-Porb re- 
lation. However, our immediate task in this section is mainly the 
construction of an updated compilation of spectroscopic SpT de- 
terminations for CV donors. We will also carry out a direct com- 
parison between this sample and a sample of isolated MS stars. 
This will allow us to establish the extent of any SpT discrepancy 
between CV donors and MS stars as a function of orbital period. 

3.1 An Updated Sample of Spectral lypes for CV 
Secondaries 

B98 presented a compilation of 54 spectroscopically established 
donor SpTs for CVs with Porb < 12 hrs. Since then, many im- 
proved and new SpT estimates have become available. We there- 
fore felt it was worth repeating their analysis in order to increase the 
reliability and size of the donor SpT sample. Briefly, we extracted 
a list of all SpT estimates for CVs with Porb < 12 hrs from the 
latest version of the Ritter & Kolb (2003) catalog (Edition 7.6). We 
inspected all of the original references for these SpTs and rejected 
all purely photometric SpT determinations, and also a few spec- 
troscopic estimates that we considered to be less compelling. In 
some cases, we also carried out additional literature searches and 
adjusted the SpTs and their errors to provide best-bet estimates 
based on all of the available evidence. Where no new information 
was available, we generally retained B98's SpT estimates, with 
two noteworthy exceptions: we were unable to find any reliable 
spectroscopic SpT determinations for OY Car and QZ Aur in the 
literature. We therefore removed these objects from our database. 
Our final compilation is provided in Table 2 and contains 91 SpTs. 

3.2 A Benchmark Main Sequence Sample 

In order to assess the departure of CV donors from thermal equi- 
librium, we need a comparison sample of SpTs for MS stars. We 
use the [A4/H] ~ sample compiled by Beuermann et al. (1999; 
hereafter B99) for this purpose. This is based on source data from 
Leggett et al. (1996; B99 Table 2) and Henry & McCarthy (1993; 
B99 Table 3). Three unresolved binaries were removed from the 
sample, and one missing spectral sub- type (GD165B: SpT=L4) 
was added from Leggett et al. (2001). In order to better cover the 
late-M and L dwarf regime, we also added a few late-type, solar 
metallicity objects from Leggett et al. (2(K)1). These again included 
some unresolved binaries, but since data is so sparse in this regime 



we did not reject these objects outright. Instead, we use these bina- 
ries only for calibrating the SpT — (J — K) relation (see below), 
where blending of objects with similar SpT should not introduce 
any serious errors. Finally, we also added the Sun (SpT=G2) and 
G134A (SpT=G3) as calibrators, in order to extend the MS sample 
into the G-dwarf regime. SpTs and absolute magnitudes for these 
two objects were taken from B98. Overall, however, we prefer the 
B99 MS sample to that used by B98 since probable low-metallicity 
objects have been removed from the former. 

In addition to providing an empirical comparison sample for 
CV donors, the MS sample is also useful for testing and calibrating 
the stellar models we will use in constructing our semi-empirical 
donor sequence. Of particular importance in this context is the abil- 
ity of the models to reproduce the observed I-K colours along the 
MS, since we will follow B98 in using this colour to estimate SpTs 
for the models. ^ 

So how well do up-to-date stellar models reproduce the I-K 
colours of MS stars? In Figure 5 we show the Mk vs (I-K) CMD 
for our MS sample, along with two sets of MS models. The first 
is the 5 Gyr isochrone taken from Baraffe et al. (1998; hereafter 
BCAH98), with colours calculated from the NextGen model atmo- 
spheres (Hauschildt, Allard & Baron 1999). The second, which we 
call AMES-MT, is an updated version of the BCAH98 models and 
has been suggested to provide a somewhat better match to observed 
optical colours (Allard, Hauschildt & Schwenke 2000). The main 
difference between the BCAH98 and AMES-MT models is the in- 
clusion of an updated TiO line list in the latter. ^ 

Somewhat surprisingly. Figure 5 shows that the older 
BCAH98 isochrone does a much better job of fitting the empiri- 
cal Mk vs (I-K) CMD. We can only conclude that the AMES-MT 
models have achieved their improvement in the optical region at 
the expense of the (near-)infrared. Given the crucial importance of 
the latter region for our purposes, we use the BCAH98/NextGen 
models throughout this paper. 

We finally revisit the SpT -{I — K) relation defined by our 
MS sample. This is shown in Figure 6. Since our sample is sub- 
stantially the same as that used by B98, their third order polyno- 
mial fit remains a good description of the data over most of its 
range. However, since our new sample now includes a mid-L dwarf 
(whereas B98's calibrating sample did not extend beyond M8), it 
makes sense to redetermine the polynomial coefficients to ensure 
that the fit also describes this very late-type regime. Figure 6 shows 
our new calibration, which is given by 

X = 56.92-35.21(7- j£:)-h9.810(7-/r)^-0.9450(7-/r)*,(13) 

where SpT = L(10 - X) for X ^ 10, SpT = M(20 - X) 
for 10 < X ^ 20, SpT = K{28 - X) for 20 < X < 28, 
and SpT = G(38 - X) for 28 < X < 38. Note that in order 
to make room for L-dwarfs, we have shifted the SpT-X relation 
(and hence the constant term of the polynomial fit) relative to that 
used by B98. As expected, our fit parameters differ only slightly 

5 Throughout this paper, we give UBV on the Johnson system, RI on the 
Cousins system, and JHK on the CIT system. 

^ Since the BCAH98 standard sequence does not provide a good match to 
solar parameters (see BCAH98), we only use these models up to 0.8 Mq 
and then supplement them witli a solar-calibrated 1.0 Mq model (taken 
from Table 3 in BCAH98). As discussed by BCAH98, the solar calibration 
is achieved by slightly varying the Helium abundance and mixing length 
parameter. Models at lower masses are not sensitive to these parameter 
changes. Note that both the AMES-MT sequence and our semi-empirical 
CV donor sequence terminate at O.TM©. 
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Table 2. Spectral types of CV donors. This table is an update of the compilation provided by B98. It is based, in the first 
instance, on the SpT estimates contained in the latest version (Edition 7.6) of the Ritter & Kolb (2003) catalog. However, 
all original source material was inspected, and only spectroscopicaUy determined SpTs that were deemed reliable were 
accepted. In some cases (e.g. where more than one SpT estimate was available in the literature), the adopted SpTs and/or 
errors were adjusted to provide best-bet estimates based on all of the available evidence. 



System P^rb (hrs) Spectral Type System Pgrf, (hrs) Spectral Type 



RXJ1951 


11.808 


MO ± 0.5 


GYCnc 


4.211 


M3 ± 1 


UYPup 


11.502 


K4±2 


SDS J2048 


4.200 


M3 ± 1 


V442 Cen 


11.040 


G6±2 


VI 043 Cen 


4.190 


M2.5 ± 0.5 


DXAnd 


10.572 


K0± 1 


SDSS J0924 


4.056 


M3.5 it 1 


AEAqr 


9.880 


K4± 1 


DO Dra 


3.969 


M4.25± 0.7 


IRXS J1548 


9.864 


K2±2 


UUAql 


3.925 


M4it 1 


ATAra 


9.012 


K2 ± 0.5 


CNOri 


3.917 


M4it 1 


RUPeg 


8.990 


K2.5 ± 0.5 


KTPer 


3.905 


M3.3 ± 1 


GYHya 


8.336 


K4.5 ± 0.5 


CYLyr 


3.818 


M3.25± 1.25 


CHUMa 


8.236 


K6.5 ± 1.5 


VYFor 


3.806 


M4.5 ± 1 


MU Cen 


8.208 


K4± 1 


IP Peg 


3.797 


M4 ± 0.5 


BT Mon 


8.012 


G8±2 


QQVul 


3.708 


M4 ± 0.5 


V1309 Ori 


7.983 


M0.5 it 0.5 


WYSge 


3.687 


M4it 1 


V392 Hya 


7.799 


K5.5 ± 0.5 


RXJ0944 


3.581 


M2it 1 


AFCam 


7.776 


K5.5±2 


MNHya 


3.390 


M3.5 it 0.5 


V363 Aur 


7.710 


G7±2 


V1432 Aql 


3.366 


M4 it 0.5 


RYSer 


7.222 


K5± 1 


TT Ari 


3.196 


M3.5 ± 0.5 


ACCnc 


7.211 


K2± 1 


MVLyr 


3.190 


M5 ± 0.5 


EM Cyg 


6.982 


K3 ± 1 


SDSS J0837 


3.180 


M5 ± 1 


Z Cam 


6.956 


K7±2 


AM Her 


3.094 


M4.25± 0.5 


SDSS J0813 


6.936 


K5.5 ± 1 


WXLMi 


2.782 


M4.5 ± 2 


V426 Oph 


6.848 


K5± 1 


RX J1554 


2.531 


M4it 1 


SSCyg 


6.603 


K4.5 ± 0.5 


SDSS J1702 


2.402 


M1.5it 1.1 


CW1045 


6.511 


K6.5 ± 1.5 


QSTel 


2.332 


M4.5 it 0.5 


CMPhe 


6.454 


M3.5 ± 1.5 


UWPic 


2.224 


M4.5 it 1 


TTCrt 


6.440 


K5 ± 0.75 


UZFor 


2.109 


M4.5 it 0.5 


BVPup 


6.353 


K3±2 


HUAqr 


2.084 


M4.25it 0.7 


AH Her 


6.195 


K7± 1 


DVUMa 


2.063 


M4.5 ± 0.5 


XY Ari 


6.065 


MO ± 0.5 


QZ Ser 


1.996 


K4 it 2 


LLLyr 


5.978 


M2.5 ± 1.5 


ARUMa 


1.932 


M5.5 ± 0.5 


AHEri 


5.738 


M4± 1 


STLMi 


1.898 


M5.5 it 1.5 


V347 Pup 


5.566 


M0.5 ± 0.5 


MR Ser 


1.891 


M6it 1 


RWTri 


5.565 


M0± 1 


CUVel 


1.884 


M5it 1 


EZDel 


5.362 


M1.5 ± 0.5 


V2301 Oph 


1.883 


M5.5 it 1 


CZOri 


5.254 


M2.5 ± 1.5 


RZLeo 


1.836 


M5it 1 


ARCnc 


5.150 


M5±l 


ZCha 


1.788 


M5.5 it 0.5 


EX Dra 


5.038 


M1.5 ± 0.5 


HTCas 


1.768 


M5.4 ± 0.3 


RXAnd 


5.037 


K4.75± 2 


V834 Cen 


1.692 


M6.5 ± 1.5 


ATCnc 


4.826 


K7.5 ± 1 


VVPup 


1.674 


M6.5 ± 1 


DQHer 


4.647 


M3 ± 0.5 


EX Hya 


1.638 


M4± 1 


AITri 


4.602 


M2.5 ± 1 


HSCam 


1.637 


M5it2 


Leo 7 


4.483 


M3±l 


BZUMa 


1.632 


M5.5 it 0.5 


SDSS J1553 


4.391 


M4.5 ± 1 


VYAqr 


1.514 


M9.5 it 1 


SS Aur 


4.387 


M3±l 


BCUMa 


1.512 


M6.5 it 0.5 


TWVu- 


4.384 


M5.5 ± 0.5 


EIPsc 


1.070 


K5it 1 


UGem 


4.246 


M4.25± 0.5 









period on the basis of its mass and radius via the period-density 
relation for Roche-lobe-filling stars; thus the assigned Port is the 
period of a hypothetical semi-detached system in which the MS star 
is the mass donor. 

This calculation requires masses and radii for all of the stars 
in our MS sample. We estimate masses from the theoretical K-band 
mass-luminosity relation - M{Mk) - predicted by the BCAH98 
models. In principle, it would be preferable to use an empirical cal- 
ibration for this purpose, in order to avoid artificially forcing the 
data points close to the models in the Porb — SpT plane. How- 



from those derived by B98, but the new fit does indeed provide an 
improved match to the latest calibrators. The RMS scatter around 
the polynomial fit is less than 0.5 sub-types for LO < SpTs < K5 
but increases to 1-2 sub-types outside this range. 



3.3 CV Donors vs Main Sequence Stars 

In Figure 7 we compare the SpTs of CV donors and MS stars. We 
follow B98 in presenting this comparison in the SpT-Porb plane. 
For this purpose, each star in the MS sample is assigned an orbital 
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Figure 5. The I-K colour -magnitude diagram for the solar-metaUicity main 
sequence sample. The solid Une shows the 5 Gyr BCAH98 isochrone, which 
is based on the NextGen atmosphere models of Hauschildt, AUard & Baron 
(1999). The dashed Une corresponds to a 5 Gyr isochrone based on a later 
version of these models (AMES-MT), which makes use of an updated TiO 
line list (Allard, Hauschildt & Schwenke 2()0()). Note that the BCAH98 
models actually provide a much better match to the data. 



ever, the best empirical M{Mk) relation (Delfosse et al. 2000) is 
only calibrated over the range 4.5 < Mr < 9.5, which excludes 
the earUest and latest stars in our sample. Moreover, Figure 3 in 
Delfosse et al. (2000) shows that the theoretical and empirical re- 
lations agree very closely in the well-calibrated regime, but there 
are hints that beyond this range, their 5th-order polynomial fit be- 
comes less reliable than the theoretical models. In any case, we 
have checked that adopting the Delfosse et al. M{Mk) relation 
(even beyond its region of validity) would not change any of our 
conclusions. 

Several radius estimates for the stars in our MS sample are 
available from B99 and Leggett et al. (2001). The bulk of our sam- 
ple comes from B99, who provide three radius estimates for each 
star. These are based on (i) the Sk [V — K) relation for MS stars, 
where Sk is the K-band surface brightness (see also Section 4.3); 
(ii) the Sk{I - K) relation; (iii) the R2{Mk) relation. B99 show 
that all of these estimates are probably good to better than 5%. We 
adopt the average of these estimates and use their range as a mea- 
sure of the associated error on the radius (and hence on the assigned 
orbital period). 

Figure 7 also shows the SpT-Porb relation predicted by the 
BCAH98 5 Gyr isochrone. For this purpose, periods were assigned 

to the models in the same way as for the MS stars (i.e. based on 
the model masses and radii), and SpTs, were estimated from Equa- 
tion 13. 

Several important results emerge immediately. First, the SpT- 
Porb relation of MS stars is fairly well described by the BCAH98 
model sequence. Second, the SpTs, of CV donors are systemati- 
cally later than those of isolated MS stars. Importantly, this is the 
case across the entire Porb range, not just for long-period systems 
with Porb ^ 3 hrs, as suggested by B98. (With hindsight, we think 



Figure 6. The spectral type vs (I-K) relation for the solar metallicity main 
sequence sample. The line shows the third order polynomial fit given by 
Equation 13. 

that a discrepancy between the SpTf, of CV donors and MS mod- 
els at short periods was aheady somewhat noticeable from Figs. 
4 and 5 of B98.) Third, apart from a few outliers, the CV donors 
with Porb ^5 — 6 hrs define a remarkably consistent SpT-Porb 
sequence. This is again evidence that most CVs do indeed follow 
a standard evolutionary track with relatively little scatter. Fourth, 
at longer periods, Porb ^5 — 6 hrs, the SpT scatter increases 
markedly, although the donors remain cooler than MS stars at the 
same period. 

B98 already provided a promising explanation for the appear- 
ance of the SpT-Porb diagram at Porb ^5 — 6 hrs. More specifi- 
cally, they showed that donors that are already somewhat nuclear- 
evolved at the start of mass transfer move through just the region of 
the SpT-Porb diagram that is occupied by observed systems. The 
increased scatter in SpTs at the longest periods could then be un- 
derstood as reflecting the range of central Hydrogen abundances in 
these donors at the start of mass transfer. The viability of this idea 
was confirmed by Baraffe & Kolb (2000) and Podsiadlowski, Han 
& Rappaport (2003). In the latter work, the authors also found that 
the dominance of evolved systems above Porb ^ 5 hrs (and the 
dominance of unevolved systems at shorter periods) is actually ex- 
pected, even if evolved systems comprise only a small percentage 
of the total CV population. For the record, the semi-empirical donor 
sequence constructed below is intended to describe the unique evo- 
lution track followed by unevolved donors, and we therefore limit 
it to Porb i$ 6 hrs. 

We finally comment briefly on the three most striking outliers 
in Figure 7, namely the three systems that lie significantly above 
the standard MS track with Porb < 3 hrs. In order of decreas- 
ing period, these systems (and their spectral type references) are 
SDSS J1702-I-3229 (Szkody et al. 2004), QZ Ser (Thorstensen et al. 
2002a) and EI Psc (Mennickent et al. 2004; see also Thorstensen 
et al. 2002b). This part of the SpT-Porb diagram should only be 
populated by low-metallicity systems (B98) or again by systems 
in which the donor was already significantly evolved at the onset 
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Figure 7. The empirical spectral type - orbital period relation for CVs and main sequence stars. Blue points are CV donors from Table 2, red points are MS stars. 
The dashed line is the relation predicted by the 5 Gyr BCAH98 isochrone, with the 5pT(I-K) calibration given by Equation 13. 



of mass-transfer. At least for QZ Ser and EI Psc, the latter is the 
preferred explanation and can also account for the fact that the or- 
bital period of EI Psc {Porb = 64 min) is well below the mini- 
mum period for "normal" CVs (Thorstensen et al. 2002ab). SDSS 
J1702-I-3229 is located within the period gap (Porb ~ 2.4 hrs), and 
the SpT estimate for it is based on the TiO-band strength measured 
by Szkody et al. (2004). This system is only about 2 a above the 
"standard" CV donor sequence, so additional observations will be 
needed to confirm if this system is genuinely abnormal. 



4 A COMPLETE, SEMI-EMPIRICAL DONOR 
SEQUENCE FOR CVS 

4.1 Constructing the Donor Sequence 

We will now use our empirical M2 — R2 and SpT-P^rb relations to 
construct and validate a complete, semi-empirical donor sequence 
for CVs. More specifically, the idea is to construct the sequence by 
combining the M2 — R2 relation with MS isochrones and atmo- 
sphere models and to validate it by checking its ability to match the 
observed SpT-Porb relation. 

The final ingredient that is needed to carry out this program 
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is a relation between M2 and T^ff (and tience SpT) along the 
donor sequence. An obvious first guess would be to assume that 
CV donors continue to follow the mass-lurrunosity relationship de- 
fined by MS stars (c.f. Patterson et al. 2003). However, this turns 
out to be incorrect. Instead, Kolb, King & Baraffe (2001; see also 
Stehle, Ritter & Kolb 1996; Baraffe & Kolb 2000; Kolb & Baraffe 
2000) show that unevolved, solar-metallicity donors should be ex- 
pected to have the same effective temperature as MS stars of iden- 
tical mass, with essentially no dependence on the mass-loss rate. 
7 

It is now straightforward to construct our semi-empirical 
donor sequence. Starting from the empirical M2 — R2 relation, we 
can use the period-density relation for CV donors to obtain Port for 
any given donor mass. The corresponding effective temperature can 
be found by interpolating on the mass-Te/ / relation in a standard 
MS isochrone. The surface gravity is of course also known (from 
M2 and R2), so absolute magnitudes in any photometric band can 
be obtained by interpolating on T^ff and log <; in a grid of model 
atmospheres and scaling to R2 . Finally, the SpT of each donor can 
be estimated from the SpT (I — K) calibration (Equation 13). 

In practice, we use the BCAH98 5 Gyr isochrone and the 
corresponding NextGen model atmosphere grid down to Te/ / ~ 
2000 K. At even lower temperatures (in the brown dwarf regime), 
the treatment of dust in the atmospheres becomes important. Be- 
tween 1500 K < Te// < 2000 K we use the 1 Gyr "DUSTY" 
isochrone of Baraffe et al. (2002; see also Chabrier et al. 2000) and 
the "AMES-DUSTY" atmosphere models of AUard et al. (2001). 
At even lower temperatures, we use the 1 Gyr COND isochrone of 
Baraffe et al. (2003) and the AMES-COND atmosphere models of 
AUard et al. (2001). Physically, these two sets of models differ in 
that dust is assimied to be present in DUSTY atmospheres, but is 
assumed to have condensed out and settled in COND models. The 
nature and location of the switch from DUSTY to COND condi- 
tions is poorly understood even in isolated brown dwarfs, and rep- 
resents a significant element of uncertainty for our donor models. 
We use younger models to represent brown dwarf donors since nu- 
clear processes only stop in CV secondaries once they have reached 
M2 ~ O.O7M0. They are therefore effectively "bom" as brown 
dwarfs at this point. The particular choice of 1 Gyr isochrones to 
represent our brown dwarf donors is based on a comparison of the 
M2 — Teff relation to the CV donor models of Kolb & Baraffe 
(1999). Isochrones at 1 Gyr do a reasonable job at reproducing this 
relation down to the lowest masses. 

We start our semi-empirical donor sequence at M2 < 
O.7OM0 (or equivalently Porb ^ 6.0 hrs). As explained in Sec- 
tion 3.3, the longer-period CV population is likely to be dominated 
by systems with evolved donors, and the evolution tracks of such 
systems depend on the central Hydrogen abundance at the onset 
of mass-transfer. By contrast, the donor sequence constructed here 
is meant to describe the unique evolution track appropriate to un- 
evolved secondaries, as are likely found in the vast majority of CVs. 
This particular starting location is also convenient since at this point 
the donor mass-radius relation just crosses the MS (Figure 3).* 

The complete donor sequence constructed in this way, span- 
ning the range O.OIM© ^ M2 < O.TOM© is hsted in Table 3 and 

^ This statement is valid down to at least M2 = 0. IMq . Like most aspects 
of the donor physics, it becomes less robust as we approach the brown dwarf 
regime (Baraffe & Kolb 2000). 

^ This cross-over is actually the correct behaviour for unevolved donors, 
since more massive stars are predominantly radiative and thus contract in 
response to mass-loss. 



plotted as a function of Port in Figure 8. More specifically, the fig- 
ure shows the evolution of the donor's physical parameters and also 
of its absolute optical and near-infrared magnitudes. For reference, 
we also show in Figure 8 the expected absolute magnitude of the 
accrction-hcatcd white dwarf as a function of Port- The effective 
temperature of the WD has been calculated following the prescrip- 
tion of Townsley & Bildsten (2003), assuming a standard theoreti- 
cal CV evolution track (Rappaport, Verbunt & Joss 1983) to yield 
M — Porb- The corresponding absolute magnitudes were obtained 
by interpolating on the WD models of Bergeron et al. (1995). These 
WD tracks provide a robust lower limit on the accretion light, since 
in reality the accretion disk is likely to dominate over the WD in 
the optical and infrared regions. 

One of the interesting aspects of the semi-empirical donor se- 
quence is the fairly sharp decline in Tg/ /, L2, and optical/infrared 
brightness in the short-period regime, but still prior to systems 
reaching Pmin- Beyond Pmin, this trend of course accelerates even 
more. In our donor track, even the transition point between L and 
T spectral types occurs before period bounce, with spectral type L 
being found only in a narrow range of periods (1.3 hrs Porb J; 
1.4 hrs). It should, of course, be kept in mind that our SpT{I-K) 
calibration is relatively poorly constrained in the L/T-dwarf regime, 
that the MS-based M2-Teff calibration becomes more inaccurate 
at the lowest masses, and that the atmosphere models become 
increasingly unreliable at the coolest temperatures. For example, 
with our caUbration, L-dwarfs are found at 1800 K < Te/ / < 
2200 K, a somewhat narrower range than that typically quoted 
for isolated L-dwarfs (1500 K < Te// < 2500 K; Kirkpatrick 
2005). On the other hand, there are also significant differences be- 
tween the effective temperatures derived for isolated L-dwarfs by 
different methods. Thus Leggett et al. (2001) foimd an L-dwarf 
temperature range of 1800 K < Te/ / < 2000 K from spectral fit- 
ting, but 1500 K < Te// < 2200 K based on the estimated radii 
and luminosities for their sample. 

Figure 8 also shows that, in period bouncers, even the 
accretion-heated white dwarf alone outshines the donor in all op- 
tical bands and is comparable to the donor in the infrared. All of 
this may explain why attempts to detect brown dwarfs in CVs (and 
especially those in suspected period boimcers) have met with very 
limited success to date. 

4.2 Validating tlie Donor Sequence 

So how well does our donor sequence fit the observed SpT-Porb 
data for CVs? The answer is provided graphically in Figure 9. In 
our view, the match to the data is very good. The average offset 
between the donor sequence and the filled data points in Figure 9 is 
consistent with zero (0.1 spectral sub-types), and the RMS scatter 
of the data around the sequence is 0.9 sub-type. Most of this scatter 
can be accounted for by observational uncertainties on the SpT 
estimates: the mean (median) SpT errors amongst the filled data 
points is 0.86 (1) spectral sub-types. For comparison, the average 
offset of the data from the standard MS is 1.2 spectral sub-types. 

It is worth stressing that the predicted sequence is not a fit 
to the data, i.e. the good match between predicted and observed 
SpTs has been achieved without any adjustable parameters. What 
this means is that the empirical M2 — R2 relation yields just the 
right amount of radius expansion to account for the late spectral 
types of the donors. Thus the results of P05 and B98 are mutually 
consistent. 

We also show in Figure 9 the donor sequence that would re- 
sult if we had adopted the MS mass-luminosity relation for CV 




Figure 8. Physical and photometric parameters along the CV donor sequence. In the left column, we show the physical donor parameters (M2, T^ff^2^ ^2 -L2) a function of 
orbital period. In the middle coltmin, the solid Unes show the optical absolute magnitudes (Mjy, Mb, My and Mji) as a function of Porb- Finally, in the right column, the sohd lines 
correspond to the red optical and near-infrared absolute magnitudes (Mj , Mj , Mh and M^- ). The wiggles in the photometric parameters near period minimum are due to switches in the 
model atmosphere grids used (BCAH98, DUSTY, COND). In the middle and right columns, we also showestimates of the absolute magnitude sequences expected for the accretion-heated 
white dwarfs in CVs (dashed Unes); see text for details. 
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Table 3. The semi-empirical donor sequence for CVs. Orbital periods are in hours, masses, radii and luminosities in solar units, and effective temperatures 
in Kelvin. UBVRI magnitudes are given on the Johnson-Cousins system (BesseU 1990), JHK magnitudes are on the CIT system (Ellas et al. 1982ab). The 
sequence provided here is abbreviated. A more complete sequence, using steps of O.OOIM© and including the far-infrared L, L' and M bands is available in 
electronic form. 



^orh 


M2 


R2 




logs 


log L2 


SpT 


Mu 


Mb 


My 


Mr 


Mj 


Mj 


Mh 


Mk 


1.462 


0.030 


0.095 


1009 


4.964 


28.50 


T 


37.82 


31.40 


26.42 


21.89 


19.00 


15.14 


15.19 


15.17 


1.420 


0.035 


0.098 


1140 


5.003 


28.74 


T 


35.32 


29.55 


25.76 


21.33 


18.37 


14.42 


14.51 


14.41 


1.384 


0.040 


0.100 


1271 


5.037 


28.96 


T 


33.52 


28.09 


25.02 


20.75 


17.80 


13.89 


13.89 


13.73 


1.353 


0.045 


0.103 


1407 


5.067 


29.15 


T 


32.22 


26.92 


24.25 


20.09 


17.13 


13.38 


13.32 


13.20 


1.326 


0.050 


0.105 


1543 


5.094 


29.33 


T 


31.13 


25.89 


23.39 


19.32 


16.42 


12.91 


12.79 


12.72 


1.302 


0.055 


0.107 


1609 


5.118 


29.42 


T 


30.83 


25.64 


23.03 


19.09 


16.27 


12.78 


12.53 


12.41 


1.281 


0.060 


0.109 


1675 


5.141 


29.51 


T 


31.68 


26.67 


23.09 


19.71 


17.21 


13.27 


12.27 


11.61 


1.268 


0.063 


0.110 


1732 


5.153 


29.57 


T 


31.22 


26.41 


22.63 


19.58 


17.21 


13.21 


12.04 


11.24 


1.290 


0.065 


0.113 


1784 


5.148 


29.65 


L6.7 


29.87 


25.28 


21.79 


18.91 


16.56 


12.74 


11.76 


11.10 


1.334 


0.070 


0.118 


1894 


5.139 


29.79 


L2.6 


28.90 


24.42 


21.15 


18.36 


15.97 


12.32 


11.46 


10.90 


1.377 


0.075 


0.123 


2003 


5.131 


29.93 


L3.7 


28.98 


23.44 


21.26 


18.68 


15.86 


11.62 


10.90 


10.66 


1.418 


0.080 


0.128 


2313 


5.124 


30.21 


M9 


24.94 


20.92 


18.94 


17.02 


14.46 


10.93 


10.26 


9.96 


1.457 


0.085 


0.133 


2477 


5.116 


30.36 


M7.6 


22.86 


19.55 


17.72 


16.05 


13.69 


10.59 


9.95 


9.63 


1.496 


0.090 


0.138 


2641 


5.110 


30.51 


M6.8 


21.03 


18.29 


16.59 


15.10 


12.96 


10.28 


9.66 


9.34 


1.533 


0.095 


0.143 


2726 


5.103 


30.59 


M6.5 


20.11 


17.64 


16.00 


14.60 


12.57 


10.10 


9.48 


9.17 


1.569 


0.100 


0.148 


2812 


5.097 


30.67 


M6.3 


19.29 


17.04 


15.44 


14.11 


12.20 


9.93 


9.32 


9.00 


1.638 


0.110 


0.157 


2921 


5.086 


30.79 


M5.9 


18.30 


16.29 


14.75 


13.50 


11.72 


9.67 


9.07 


8.76 


1.704 


0.120 


0.166 


2988 


5.076 


30.88 


M5.6 


17.70 


15.81 


14.30 


13.10 


11.40 


9.48 


8.88 


8.58 


1.768 


0.130 


0.175 


3055 


5.066 


30.96 


M5.3 


17.16 


15.38 


13.90 


12.74 


11.12 


9.29 


8.71 


8.40 


1.828 


0.140 


0.183 


3103 


5.057 


31.03 


M5 


16.76 


15.06 


13.59 


12.46 


10.89 


9.14 


8.55 


8.26 


1.886 


0.150 


0.192 


3151 


5.049 


31.10 


M4.8 


16.37 


14.74 


13.30 


12.20 


10.68 


9.00 


8.42 


8.12 


1.942 


0.160 


0.200 


3185 


5.042 


31.15 


M4.6 


16.08 


14.50 


13.07 


11.99 


10.51 


8.87 


8.30 


8.00 


1.996 


0.170 


0.207 


3218 


5.035 


31.20 


M4.4 


15.82 


14.28 


12.86 


11.80 


10.36 


8.75 


8.18 


7.89 


2.049 


0.180 


0.215 


3246 


5.028 


31.25 


M4.3 


15.60 


14.09 


12.68 


11.64 


10.22 


8.65 


8.07 


7.79 


2.100 


0.190 


0.223 


3269 


5.021 


31.29 


M4.2 


15.42 


13.93 


12.53 


11.50 


10.10 


8.55 


7.98 


7.69 


2.149 


0.200 


0.230 


3292 


5.015 


31.33 


M4 


15.24 


13.77 


12.38 


11.36 


9.98 


8.45 


7.88 


7.60 


3.179 


0.200 


0.299 


3292 


4.789 


31.56 


M4.2 


14.59 


13.19 


11.85 


10.83 


9.44 


7.88 


7.31 


7.02 


3.559 


0.250 


0.347 


3376 


4.756 


31.73 


M3.8 


13.95 


12.58 


11.27 


10.29 


8.95 


7.47 


6.90 


6.63 


3.902 


0.300 


0.392 


3436 


4.729 


31.87 


M3.5 


13.48 


12.13 


10.84 


9.88 


8.58 


7.15 


6.57 


6.31 


4.218 


0.350 


0.434 


3475 


4.706 


31.98 


M3.3 


13.13 


11.79 


10.51 


9.57 


8.29 


6.89 


6.31 


6.05 


4.512 


0.400 


0.475 


3522 


4.686 


32.08 


M3.1 


12.80 


11.46 


10.19 


9.26 


8.02 


6.65 


6.06 


5.82 


4.788 


0.450 


0.514 


3582 


4.669 


32.18 


M2.8 


12.47 


11.13 


9.86 


8.96 


7.75 


6.42 


5.83 


5.60 


5.050 


0.500 


0.552 


3649 


4.653 


32.27 


M2.5 


12.16 


10.81 


9.54 


8.66 


7.49 


6.20 


5.60 


5.39 


5.299 


0.550 


0.588 


3760 


4.639 


32.38 


M2 


11.79 


10.40 


9.14 


8.28 


7.19 


5.97 


5.35 


5.16 


5.537 


0.600 


0.623 


3893 


4.626 


32.49 


M1.4 


11.42 


9.97 


8.70 


7.88 


6.88 


5.73 


5.08 


4.93 


5.766 


0.650 


0.658 


4065 


4.615 


32.61 


M0.4 


10.96 


9.45 


8.19 


7.40 


6.54 


5.48 


4.83 


4.71 


5.985 


0.700 


0.691 


4239 


4.604 


32.73 


K7 


10.46 


8.94 


7.70 


6.94 


6.20 


5.24 


4.62 


4.52 


6.198 


0.750 


0.724 


4460 


4.593 


32.85 


K5 


9.74 


8.34 


7.15 


6.43 


5.85 


5.00 


4.42 


4.35 



donors (rather than the correct mass-effective temperature one). As 
expected, this predicts even later spectral types, since in this case 
the bloated donor must reduce its Te/ / in order to maintain the MS 
luminosity that would be appropriate for its mass. This sequence 
does not match the data nearly as well - the average offset from the 
data is 1 .4 spectral sub-types - but it does provides a useful lower 
limit on the predicted SpTs. 

We finally comment again on the remarkably small scatter in 
the SpT-Porb diagram for unevolved CVs. In principle, we could 
try to estimate the intrinsic dispersion demanded by the data in the 
same way as in Section 2.3, i.e. by demanding that = 1 for 
the correct aint,SpT. This would yield aint,SpT ^ 0.4. However, 
this value should perhaps not be taken at face value, since most 
spectral type errors are ad hoc estimates that cannot be expected to 
be Gaussian or even strictly l-tr. A more robust statement is that 
the average intrinsic dispersion must be less than the RMS scatter 
of about 1 spectral sub-type. Tliis is still quite impressive: in the M- 
dwarf regime, one sub-type typically spans an effective temperature 
range of only ATe/ / ~ 100 - 250 K. 



4.3 Applying the Donor Sequence: Distances from 
Photometric Parallaxes 

We conclude our study by exploring one obvious practical appli- 
cation of the semi-empirical donor sequence, namely distance esti- 
mation. Distances towards CVs are of fundamental importance to 
virtually all areas of CV research, but are notoriously difficult to de- 
termine. In recent years, trigonometric parallaxes have finally been 
determined for a number of systems, but indirect methods remain 
the only way to obtain distance estimates for the vast majority of 
CVs. 

Perhaps the most widely used indirect technique for estimat- 
ing distances towards CV is Bailey's (1981) method. This is based 
on the K-band surface brightness (Sk) of the secondary, defined as 

=if + 5-51ogd + 51ogi?2/i?0 (14) 

where d is the distance in parsecs. Sk can be calibrated against 
colour (usually {V — K)) by using a suitable sample of MS stars. 
The distance towards a CV can then be estimated if the apparent K- 
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Figure 9. Observed vs predicted spectral types of CV donors. Points correspond to empirically determined SpTs for CVs with Porb ^ 6 hrs, taken 
from Table 2. Likely evolved donors are shown as open symbols and were excluded from the statistics quoted in the text. The solid red line shows the 
SpT(Porb) relation predicted by the semi-empirical donor sequence. The dashed line shows the relation predicted by the BCAH98 5 Gyr isochrone. 
The dotted line shows the relation that would be obtained if CV donors obeyed the MS mass-luminosity relation, rather than the MS mass-effective 
temperature relation. The kinks in the solid and dotted lines are due to a switch between two different model atmosphere grids. 



band magnitude, V — K colour and radius of the donor are known. 
In practice, the K-band magnitude of the CV is usually assumed 
to be dominated by the donor (so the resulting distance estimate is 
really a lower limit). Also, the radius is generally estimated from 
the orbital period, by using the Porb-P2 relation (Equation 3) and 
assuming a MS mass-radius relation. However, the most difficult 
aspect of the method is the determination of an appropriate V — K 
value for the secondary, since the optical flux of CVs is usually 
dominated by accretion light. 

Given this difficulty, much of the original promise of Bailey's 
method derived from the form he inferred for the Sk{V — K) re- 



lation. He found that for M-dwarfs (and thus for almost all CV 
donors), Sk was essentially constant, so that even large errors in 
V — K would not seriously affect the resulting distance estimates. 
However, the ^^(V^ — K) relation has since been recalibrated by 
Ramseyer (1994), B99 and, most importantly, Beuermann (2000). 
The last of these studies provides by far the cleanest calibration to 
date and shows definitively that Sk is not constant in the M-dwarf 
regime. Instead, the Sk{V — K) relation is approximately linear 
over essentially the full range of SpTs found in CVs. As a result. 
Bailey's original method is much less robust than had previously 
been supposed, unless the contribution of the donor to both V and 
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K band fluxes can be estimated reliably from observations. This is 
usually impossible, and the only way forward is then to assume a 
"typical" donor {V — K) for any given CV, based perhaps on its 
observed SpT or just on its orbital period. 

Our semi-empirical donor sequence provides a way to sim- 
plify and improve such distance estimates. As already noted above, 
the small scatter around the M2 — R2 and the SpT — Port rela- 
tions implies that most CVs do, in fact, follow the unique evolution 
track that is delineated by this sequence. We can therefore use the 
absolute magnitudes predicted along the sequence to obtain lower 
limits on the distance towards any CV, under the single assump- 
tion that the system follows the standard track. For example, the 
lower limit on the distance associated with a single epoch A'-band 
measurement is 

K ~MK,2{Porb)+5 



log di, 



(15) 



where K is the apparent magnitude and Mk,2 is the absolute K- 
band magnitude on the donor sequence at the CVs orbital pe- 
riod. In principle, the apparent magnitude should be extinction- 
corrected, but in practice this correction is usually negligible for 
CVs in the infrared. If the donor contribution to the total flux is 
known, the apparent magnitude of the system should of course be 
replaced with that of the donor to yield an actual estimate of the 
distance (rather than just a lower limit). 

In order to test this method, we have compiled distances and 
2MASS infrared magnitudes for all CV with trigonometric paral- 
lax measurements. The sample is listed in Table 4 and contains 23 
systems, 22 of which have reliable 2MASS measurements. In Fig- 
ure 10, the absolute IHK magnitudes are shown as a function of 
orbital period and compared to the predictions for the donor stars 
from the semi-empirical sequence in Table 3. We have intentionally 
not corrected the observed data points for interstellar extinction. In 
practical applications, extinction estimates will often not be avail- 
able, and our main goal here is to test how well distances may be 
estimated just from single-epoch infrared magnitudes and orbital 
periods. Extinction estimates are nevertheless listed in Table 4 and 
are negligible for the sources in our sample. ^ 

Perhaps the most important point to note from Figure 10 is 
that the semi-empirical donor sequence does a nice job of tracing 
the lower envelope of the data points. No CV with Porb < 6 hrs is 
significantly fainter in the infrared than the donor prediction. This 
again validates the sequence and implies that it can be used with 
confidence to determine lower limits on CV distances. However, 
it is also obvious from Figure 10 that the majority of systems are 
significantly brighter than the pure donor prediction, even in the 
K-band. Thus in most CVs, even the infrared flux is dominated by 
accretion light (cf. Dhillon et al. 2000). 

We have determined the average offsets between the donor se- 
quence and the data in the three infrared bands, for systems with 
Porb < 6 hrs. In doing so, we assumed that none of the systems 
in our sample are period bouncers, except WZ Sge, which was re- 
moved from the sample for the purpose of this calculation. The 
resulting offsets are 



AJ = 1.56 
AH = 1.34 
AK = 1.21 



(TAJ = 1.25, 
(TAJ = 1.12, 
(TAJ = 1.07, 



(16) 



^ There is of course a significant bias here, in that our trigonometric par- 
allax sample is, by definition, nearby and thus less heavily reddened than a 
more "typical" CV sample. However, infrared extinction is very unlikely to 
ever dominate the error budget, even for more distant CVs. 



J-Band 



H-Band 




K-Band 



Period (hrs) 



Figure 10. Infrared absolute magnitudes of CVs as a function of Porb- 
Points con'espond to absolute magnitudes for CVs with trigonometric par- 
allaxes and reliable 2MASS observations (Table 4). Colours correspond to 
various CV sub-types, as indicated (DN = dwarf nova; SU = SU UMa star; 
NL = non-magnetic nova-like; AM = AM Her star; IP = intermediate po- 
lar). Systems with Porb > 6 hi's are shown as open symbols. Solid lines 
show the tracks predicted by the semi-empirical donor sequence, and dotted 
lines show the tracks corresponding to the BCAH98 5 Gyr isochrone. The 
dashed lines show the CV donor tracks offset by the average offset between 
the data points and the pure donor predictions; see text for details. 



where ctajhk is the scatter of the data around the mean offsets. 
As one might expect, both the offsets and the scatter decrease with 
increasing wavelength. However, even at K, the donor typically 
contributes only 33% of the total flux. The corresponding contri- 
butions at J and H are 24% and 29%, respectively. Thus a distance 
limit obtained from a single epoch infrared measurement with no 
additional information would typically underestimate the true dis- 
tance by factors of 2.05 (J), 1.86 (H) and 1.75 (K). 

If it was deemed important to convert these robust lower lim- 
its into (much less robust) distance estimates, one might consider 
applying the mean offsets to the absolute donor magnitudes. We 
hasten to emphasize that this procedure has no physical basis, and 
indeed there is no reason to think a priori that a single offset should 
yield reasonable results across the full range of orbital periods 
and CV types. Nevertheless, the effect of applying these offsets to 
the donor sequence is illustrated in Figure 10 (dashed lines). The 
o'jHK values indicate that, for our sample, the resulting distance 
estimates would be correct to about a factor of 1.78 (J), 1.68 (H) 
and 1.64 (K). 
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Table 4. Distance moduli and absolute infrared magnitudes for CVs with trigonometric parallaxes. Apparent magnitudes were 
taken from the final 2MASS catalog and converted to the CIT system using the transformations given by Carpenter (2001). Note 
that HV Vir has a parallax but no reliable 2MASS detection. Note that any asymmetric errors on the distance moduli have been 
symmetrized for simplicity; see the original references for the exact imcertainties. We also provide estimated reddening values 

(Eb-v) from Verbunt (1987). For reference, the associated extinctions are Aj = 0.89 X Eb-v, Ah = 0.43 X Eb-v and 
Ak = 0.,35 X Eb-v (Rieke & Lebofsky 1985). All systems that are not included in Verbunt (1987) are very close and are 
therefore listed here with reddening entries of (0.0). The distance references indicated in the last column are 1 = Duerbeck (1999); 
2 = Harrison et al. (2004); 3 = Thorstensen (2003); 4 = McArthur et al. (1999); 5 = McArthur et al. (2001); 6 = Beuermann et al. 
(2004); 7 = Beuermann et al. (2003). 



System 


Porb{hrs) 


Type 


m — M 


Mj 


Mh 


Mk 


Eb-v 


Ref 


AE Aqr 


9.880 


no 


5 04 + 65 






3 76 + 65 


(0.00) 


1 


RU Peg 


8.990 


DNAJG 


7.29 ± 0.17 


^ 7Q _|_ 17 
J. / o ^ yj. 1 / 


'2 -2^ _|_ 17 
J, J J ZL \J, 1 / 


'2 90 -1- n 17 


0.00 
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ZCam 


6.956 


DN/ZC 


6.06 ± 0.67 


5.51 ± 0.67 


4.98 ± 0.67 


4.82 ± 0.67 


0.00 


3 


SS Cyg 


6.603 


DN/UG 


6.13 ± 0.15 


2.41 ± 0.15 


2.22 ± 0.15 


2.19 ±0.15 


0.04 


2 


AH Her 


6.195 


DN/ZC 


9.10 ± 0.76 


2.89 ± 0.76 


2.63 ± 0.76 


2.58 ± 0.76 


0.03 


3 


RWTri 


5.565 


NL 


7.79 ± 0.22 


4.16 ±0.22 


3.78 ± 0.22 


3.69 ± 0.22 


0.10 


4 


TV Col 


5.486 


IP 


7.87 ± 0.10 


5.34 ±0.10 


4.95 ± 0.10 


4.85 ±0.10 


0.00 


5 


V3885 Sgr 


4.971 


NL 


5.21 ± 0.48 


4.76 ± 0.48 


4.52 ± 0.48 


4.43 ± 0.48 


0.00 


1 


IXVel 


4.654 


NL 


4.91 ± 0.20 


4.23 ± 0.20 


4.06 ± 0.20 


3.94 ± 0.20 


0.00 


1 


SS Aur 


4.387 


DNAJG 


6.13 ±0.13 


6.57 ±0.13 


6.13 ±0.13 


5.89 ±0.13 


0.10 


2 


UGem 


4.246 


DNAJG 


5.01 ± 0.09 


6.63 ± 0.09 


6.05 ± 0.09 


5.84 ± 0.09 


0.00 


2 


V1223 Sgr 


3.366 


IP 


8.61 ± 0.20 


4.23 ± 0.20 


4.12 ±0.20 


4.05 ± 0.20 


0.15 


6 


AM Her 


3.094 


AM 


4.49 ±0.19 


7.21 ±0.19 


6.70 ±0.19 


6.54 ±0.19 


(0.00) 


3 


YZCnc 


2.083 


DN/SU 


7.54 ± 0.29 


5.64 ± 0.29 


5.40 ± 0.29 


5.31 ± 0.29 


0.00 


2 


SUUMa 


1.832 


DN/SU 


7.07 ± 1.06 


4.73 ± 1.06 


4.65 ± 1.06 


4.62 ± 1.06 


0.00 


3 


V893 Sco 


1.823 


DN/SU 


5.95 ± 0.61 


7.28 ±0.61 


6.95 ± 0.61 


6.75 ± 0.61 


(0.00) 


3 


EXHya 


1.638 


IP 


4.05 ± 0.04 


8.23 ± 0.05 


7.89 ± 0.05 


7.66 ± 0.05 


0.00 


7 


VYAqr 


1.514 


DN/SU/WZ 


4.93 ± 0.30 


10.34 ± 0.30 


9.91 ±0.31 


9.68 ±0.31 


(0.00) 


3 


TLeo 


1.412 


DN/SU 


5.02 ± 0.26 


9.74 ± 0.26 


9.30 ± 0.26 


8.83 ± 0.26 


(0.00) 


3 


HVVir 


1.370 


DN/SU/WZ 


8.31 ± 1.37 










3 


WZSge 


1.361 


DN/SU/WZ 


3.19 ±0.14 


11.66 ±0.15 


11.35 ±0.15 


10.83 ±0.15 


0.00 


2 


EFEri 


1.350 


AM 


6.06 ± 0.77 


11.09 ±0.80 


9.60 ± 0.78 


9.33 ± 0.79 


(0.00) 


3 


GWLib 


1.280 


DN/SU/WZ 


5.09 ± 0.51 


11.10 ±0.51 


10.49 ± 0.52 


10.33 ± 0.54 


(0.00) 


3 



5 DISCUSSION AND CONCLUSIONS 

The main goal of this paper has been the construction of a com- 
plete, semi-empirical donor sequence for CVs that is based on our 
best understanding of donor properties and is consistent with all 
key observational constraints. This donor sequence is provided in 
Table 3 and is based primarily on an updated version of the mass- 
radius relationship for CV secondaries determined by P05. It also 
relies on a MS-based mass-effective temperature relation that has 
been shown to be appropriate for CV donors by theoretical work, 
and on up-to-date stellar atmosphere models. By design, the donor 
track also reproduces the observed locations of the period gap and 
the period minimum, and extends beyond this to the period bouncer 
regime. We have shown that this sequence provides an excellent 
match to the observed SpT-Porb relation for unevolved CV donors 
(i.e. Porb 6 hrs), and that it correctly traces the lower envelope 
of the MjHK — Port distributions of CVs. 

Along the way, we have revisited and updated two important 
results on CV donors by P05 and B98. Regarding the former study, 
we have carried out a full, independent analysis of the superhumper 
and eclipser data that was used by P05 to construct their empirical 
A/2 — R2 relation. We constructed and used our own q — e cali- 
bration for superhumpers, verified the validity of the constant Mi 
assumption on which the method rests, imposed the locations of 
the period gap and the period minimum as constraints on the mass- 
radius relation and used a self-consistent fitting technique to extract 
optimal parameters for the M2 — R2 relations for long-period CVs, 
short-period CVs and period-bouncers. 



Regarding the study by B98, we have updated their SpT data 

base for CV secondaries. Our new sample contains 91 CV donors 
with spectroscopically-determined SpT estimates (compared to 54 
donors in the B98 sample). Comparing the SpT-PoTb distribution 
to that of MS stars, we find that CV donors have later spectral types 
at all orbital periods. This extends the finding of B98 to the short- 
period regime {Porb ^ 3 hrs). 

We find that there is remarkably little intrinsic scatter about 
both the mean M2 — R2 and SpT — Porb relations for CVs with 
Porb < 6 hrs: probably no more than a few percent in R2 and less 
than 1 spectral sub-type. This suggests that, on long time-scales, 
most CVs do indeed follow a unique evolution track, just as is the- 
oretically expected (e.g. Kolb 1993; Stehle, Ritter & Kolb 1996). 
The large scatter in luminosity-based accretion rate estimates at 
fixed Porb (e.g. Patterson 1984) is therefore probably caused by 
fluctuations around the mean mass transfer rate on time-scales that 
are short compared to the time-scale on which the donor loses mass 
and the binary evolves (see Biining & Ritter 2004 and references 
therein). In principle, the empirical M2 — R2 relationship for CV 
donors can be inverted to infer the long-term average mass-transfer 
rate as a function of orbital period; this will be the subject of a 
separate paper. 

An important feature of our donor sequence is a sharp decline 
in effective temperature, luminosity and optical/IR brightness, well 
before the period minimum is reached. In fact, the L/T transition es- 
sentially coincides with Pmin on the evolution track (although this 
statement depends somewhat on the validity of our spectral type 
calibration in a sparsely populated regime). Beyond Pmin, even 
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the accretion-heated WD alone is expected to outshine the donor at 
optical wavelengths and it is comparably bright to the donor even in 
the infrared. All of this helps to explain why brown dwarf donors 
have proven so difficult to detect in CVs, especially in suspected 
period bouncers. 

Finally, we have looked at one obvious application of our 
donor sequence, the determination of distance estimates from in- 
frared photometry. We have shown that robust lower limits can in- 
deed be obtained, but that, in our calibration sample, these are typ- 
ically a factor of two smaller than the true distances. Thus, even 
in the infrared, the donor typically contributes only ~30% of the 
total flux. Distance estimates that explicitly allow for this average 
donor contribution are typically good to about a factor of 2 in our 
calibration sample. 
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APPENDIX A: CALIBRATING THE e - q RELATION 

Here we provide some additional details on the way in which we 
calibrated the e — g relation. As already noted in Section 2.2, we rely 
on the same calibrators as P05 (given in their Table 7), but there are 
a few differences between their treatment of the problem and ours. 
First, the two calibrators with the most precise mass ratios are XZ 
Eri and DV UMa (Feline et al. 2004). In both cases. Feline et al. 
obtained their final g-estimates as weighted averages over the mass 
ratios determined in 3 different photometric bands. However, we 
foimd the individual g-estimates for both systems were not quite 
consistent at the level expected based on their internal errors. We 
therefore re-averaged the individual g values, but allowed explic- 
itly for just enough intrinsic dispersion to make the values in all 
bands consistent with each other. The resulting new mass ratios 
were g = 0.111 ± 0.03 (XZ Eri) and g = 0.152 ± 0.003 (DV 
UMa). The implied changes to the mass ratios are relatively small, 
but the increase in the formal errors is important (since these two 
data points would otherwise be given too much weight in the cal- 
ibration). Second, we took note of the approximate upper limit on 
the mass ratio of BB Dor suggested by P05 (g < 0.38), but did not 
explicitly enforce it in our fits. Instead, we decided to check a pos- 
teriori if our fits significantly violated this approximate constraint 
or not. Third, we chose to treat g (rather than e) as the independent 
variable of the problem, since the ultimate goal is to estimate g from 
a measured e. Fourth, we experimented with a variety of different 
families of 2-parameter curves, including linear fits (with constant 
term), quadratic fits (with no constant term) and power laws. 

We found that the preferred functional form for the e — g rela- 
tion is linear, i.e. q = a + be, with a and b constants. We obtained 
optimal estimates of a and b by minimizing the of the model 
with respect to the data, allowing explicitly for errors in both g and 
e and also for the possibiUty of intrinsic scatter aroimd the relation, 
Tint - The goodness-of-fit metric then becomes 



(1- 



■bef 



(Al) 



The e — g relation given in the main body is the minimum x 
solution. It does not significantly violate the approximate upper 
limit on BB Dor's mass ratio and achieves a statistically accept- 



able xt = 103 with 



0. Thus any intrinsic scatter around 



the calibrating relation must be small compared to the statistical er- 
rors on the data point. P05 explicitly added an extra 5% uncertainty 
on both g and e when carrying out their fits. We appreciate their 
rationale for doing this (allowing for variability in superhump pe- 
riods and unaccounted-for external errors in published mass ratio 
estimates); indeed, we accoimted for the possibiUty of a non-zero 
Tint in our fits for exactly the same reasons. However, it simply 
turns out that such additional error contributions are not actually 
demanded by the data. 



APPENDIX B: TESTING THE ASSUMPTION OF 
CONSTANT PRIMARY MASS 

In order to test if there is evidence for evolution of Mi, we fitted 
a straight line to the A'^ = 16 eclipser data points with Porb < 
6 hrs given in Table 7 in P05. Orbital period errors are negligible 
compared to the errors on Mi and to the intrinsic dispersion, so we 



used the x metric 



E 



(Ml 



bPorbf 



(Bl) 



Here, ctmi is the error on Mi and CTi„t is the intrinsic dispersion 
around the fit. The intrinsic dispersion is set to the value needed 
to obtain a reduced — 1. where v = A*' — 2; for the linear 
fit, this value is <7int = O.IGM©. The corresponding constraints 
on the slope and intercept are shown in the top panel of Figure 2, 
and the best-fitting line itself is plotted in the bottom panel. There 
is certainly substantial scatter among the Mi values, but the best- 
fitting slope is not significantly different from zero. Thus there is 
no convincing evidence for a trend of Mi with Porb- The best- 
fitting constant Mi, using the same x^ metric, but with slope fixed 
at 6 = 0, yields a mean WD mass of Mi = 0.75 ± O.O5M0. 
It requires the same level of intrinsic scatter as the linear fit, i.e. 

(Tint = O.16M0. 

These results may seem at odds with the studies of Webbink 
(1990) and Smith & Dhillon (1998), both of whom suggested that 
CVs below the gap tend to have lower WD masses than those above 
the gap. More specifically. Smith & Dhillon [Webbink] derived 
weighted mean masses of Ml = O.66±O.O1M0 [O.eGiO.OlMo] 
below the gap and Mi = 0.78±0.02Mq [O.81±O.O4M0] above 
the gap. Two points are important to note here. First, the weights 
that were used in calculating these weighted means were the in- 
verse variances of the formal errors on the individual data points. 
Second, the quoted uncertainties on the means are the formal er- 
rors on the weighted means, rather than the dispersion of the points 
around the means. If we calculate the same quantities for our sam- 
ple (in which there are 8 short- and 8 long-period systems), we find 
weighted means with formal errors of Mi = 0.642 ± O.OO7M0 
below the gap and Mi = 0.72 ± O.O2M0 above the gap. So if 
we had used these measures, there would appear to be a significant 
difference between short- and long-period systems for our sample, 
too. 
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The resolution of this apparent paradox is simple. As already 
noted by Smith & Dhillon, the weighted means and their formal 
errors are quite misleading, since they can be totally dominated by 
a few systems with small formal uncertainties. If the true intrinsic 
dispersion among the sample being averaged is larger than these 
formal uncertainties - as it clearly is here - then the formal error 
on the weighted mean (and even the weighted mean itself) is rather 
meaningless. From a statistical point of view, such a sample vio- 
lates the basic premise underlying the calculation of the weighted 
mean, which is that each data point is drawn from a distribution 
with the same mean, with a variance given solely by the formal 
error. 

Smith & Dhillon also give the unweighted means and disper- 
sions for their sample, which are Mi = 0.69 ± 0.13A'/q below the 
gap and Mi = 0.80 ± O.22M0 above the gap. For our sample, the 
corresponding numbers are Mi = 0.74 ± O.2OM0 below the gap 
and Ml = 0.78 ± O.IGM© above the gap. This already shows that 
the differences between short- and long-period systems are much 
less significant than suggested by the formal weighted means and 
errors, and that any real trend with Porb that might be present is 
small compared to the intrinsic dispersion. 

However, the best way to estimate the means, error on the 
means and intrinsic dispersion is via the statistic given by Equa- 
tion Bl (again with slope fixed at fe = 0). This accounts explic- 
itly for both the formal errors and intrinsic dispersion, with the 
latter being adjusted to yield a reasonable xt. If we do this sep- 
arately for short- and long-period systems in our sample, we find 
Ml = 0.73 ± 0.07Mq with = 0.19Mq below the gap and 
Ml = 0.77 ± O.O6M0 with oint = O.ISM© above the gap. As 
expected, these numbers are not significantly different from each 
other anymore. 



APPENDIX C: FITTING THE M2 R2 RELATION 

We wish to obtain an optimal estimate of the mass-radius relation 
of CV donors from a set of (M2 , R2) pairs. As noted in Section 2.5, 
the key challenges are to acccount for the correlated nature of the 
errors on M2 and R2 and to self-consistently impose the external 
constraints derived from the locations of the period gap and the 
period minimum. In addition, we also need to allow for intrinsic 
scatter around the mass-radius relation. 

We begin by defining the analytical form of the mass-radius 
relation we wish to fit to the data. Based on inspection of the data in 
Figure 3, we will follow P05 and describe the mass-radius relation 
as a power law, 



Rrt 



R 



ref 



M 



M, 



ref 



(CI) 



Here and below, we drop the subscript "2" on the mass and radius 
in order to keep the notation transparent. In practice, we carry our 
fits out in log-space, where the power law transforms into the linear 
relation 



log Rmod = log Rraf -I- 6 log M - 6 log Mrcf ■ 



(C2) 



In the absence of additional constraints, both Rref and b 
would be free parameters of this model. However, as discussed in 
Section 2.4, we will demand that M = Mconv at Port = Pgap,± 
and M — Mbounce at Port = Pmin- The easiest way to ac- 
complish this is to adopt Mref = Mconv for long- and short- 
period CVs and Mref = Mbounce for period bouncers. With these 
choices of reference mass, Rref is fixed (to within some error) by 



the period-density relation for Roche-lobe filling stars (Equation 6). 
Thus for long period CVs, we have 



logi?,e/ = logO.2361i?0 + (2/3) logPgap,+ {hr)+ 

{l/3)logMeonv -(l/3)logMQ. 



(C3) 



Similar relations hold for short-period CVs and period bouncers. 
Note that in all these relations, Pgap,±(or Pmin) and Mconv (or 
Mbounce) are empirical estimates with associated unertainties. As 
discussed further below, these uncertainties translate into a system- 
atic error on each data point that needs to be accounted for in the 
final fit. 

However, let us first deal with the correlation between the M 
and R estimates for a given system. The easiest way to achieve 
this is to work exclusively in terms of M. For example, consider 
a particular data point (Mi, Ri, Porb,i) in the short-period regime 
(the corresponding derivation for long-period systems and period 
bouncers is identical). The residual between this point and the 
model prediction is — log Ri — log Rmod,i{Mi). Using Equa- 
tions 6, CI and C3, this can be rewritten as 



Ai = 3 log Porb,i — 

(i-6)logM, 



l0gP3ap,- + 

- (i -6) log Me, 



(C4) 



Let us ignore the systematic errors arising from the uncertainties on 



Mc. 



and Pgap,- for the moment and consider only the statistical 



variance on Ai. This is 



2 



(C5) 



where we have neglected the error on Porb,i which is always much 
smaller than that on Mi. An optimal estimate for b could then be 
obtained by minimizing the usual statistic 



., stat.i 



(C6) 



or, if we additionally want to allow for an intrinsic variance (Jint> 



(C7) 



and 



where N is the number of data points. 

We should also face up to the uncertainties on Mcon 
Pgap,-. It is tempting to simply account for these in the denom- 
inator of Equation C7. However, this would be incorrect, since any 
change in these quantities affects all residuals in the same way: 
cTiog Mco„„ and aiog Pgap,. are systematic errors in this context. Let 
us define the systematic variance on each residual as 



(3 



2 

Clog Mo, 



(C8) 



It is then possible to define a new statistic that has the cor- 
rect distribution and accounts for both (Jstat,i and CTs^s- Following 
Stump et al. (2002), we find that in our case, this statistic is 



JV 



A? 



._, stat,i ' int 



(C9) 



The only uncertainties we have not explicitly accounted for in 
our fits are (i) the error on the mean white dwarf mass adopted for 
superhumpers^" and (ii) the errors on the parameters of the £— g cal- 



Note that we are specifically referring to the error on the mean here. 
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ibration. As discussed in Sections 2.3 and 2.2, these uncertainties 
also translate into systematic errors on the data points. In principle, 
it would be possible to formally account for these systematics in a 
similar way as we just did for the errors on Pgap,± and Mconv, for 
example. However, in practice, such a treatment would be consid- 
erably more difficult and (in our view) a little pointless. The added 
difficulty arises partly because these new systematics do not affect 
all data points in a given mass/period regime (i.e. only the super- 
humpers, but not the eclipsers); the error on the calibration parame- 
ters, in particular, does not even affect all superhumpers in the same 
way. We feel the additional complexity this would introduce is not 
warranted. After all, regarding (i), we already know that the error on 
the mean Mi is small compared to intrinsic dispersion around this 
mean, which we do account for (see Section 2.3). And regarding 
(ii), we already noted in Section 2.2 that the systematic uncertainty 
associated with the e — g relation is small in the well-calibrated 
regime, but not necessarily well-described by the formal parameter 
uncertainties in the poorly-calibrated regime. Rather than attempt 
to include this poorly defined systematic in our fits, we thus prefer 
to simply emphasize its existence. 

The optimal estimate of the power-law exponent b is obtained 
by minimizing Equation C9. For long- and short-period systems, 
we determine the appropriate level of aint by requiring that the 
reduced should be equal to one; for period bouncers, the reduced 
is slightly less than one even without any intrinsic dispersion. 
However, this is probably more of a reflection of the sparseness 
of the data in this regime than of any true constraint on (7i„t- We 
therefore conservatively adopt the value of aint obtained for short- 
period CVs for the period bouncers as well. Errors on the power law 
exponents can be estimated in the usual way. For example, the 1-cr 
confidence interval around b corresponds to the range of exponents 
for which ^ Xmin + 1^ where xinin is the lowest value of x^- 

In order to test this method, we have carried out Monte Carlo 
simulations. Thus we created fake data sets with Imown slopes and 
subject to all of the errors we are trying to account for in our fits. 
In these simulations, the method did a good job in recovering the 
correct slope, and the distribution of the recovered slopes had a 
dispersion consistent with the estimated error on the slope. 
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